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ABSTRACT 

Recent analysis of the Milky Way's satellite galaxies reveals that these objects share 
a common central mass density, even though their luminosities range over five orders 
of magnitude. This observation can be understood in the context of galaxy formation 
theory by quantifying the factors which restrict the central mass density to a small 
range. One limit is set by the maximum mass that can collapse into a given region 
by the hierarchical growth of structure in the standard cold dark matter cosmology. 
Another limit comes from the natural thresholds which exist for gas to be able to cool 
and form a galaxy. The wide range of luminosities in these satellites reflect the effects 
of supernova feedback on the fraction of cooled baryons which are retained. 



1 INTRODUCTION 

The satellite galaxies of the Milky Way are providing a valu- 
able guide to our understanding of cosmology. Though they 
constitute only a single statistical sample, the standard cold 
dark matter cosmological model (ACDM) must at least be 
able to account for their abundance and properties. Attain- 
ing such consistency is not proving to be straightforward. 

Initial friction between observation and theory con- 
cerned the number of observed satellites, which is greatly 
exceeded by the number of subhalos that ACDM predicts 
in a host halo of the Milky Way's size (Klypin et al. 1999; 
Moore et al. 1999, and over 1000 citing articles). A proposed 
resolution to this was that photoionisation has suppressed 
the formation of small galaxies in these subhalos (Bullock 
et al. 2001; Benson et al. 2002; Somerville 2002). 

Observations of satellite galaxies also indicate that 
these systems boast extremely high mass to light ratios, 
much greater than are deduced for more luminous galaxies. 
Mass does not appear to follow light, as might be naively 
expected. This has been recently quantified by Strigari et 
al. (2008), leading to a particularly intriguing result. The 
central mass of each satellite, as contained within a certain 
chosen radius (300 pc), appears to be the same for all the 
satellites, remaining approximately constant over seven or- 
ders of magnitude in luminosity. 

A suggested explanation for this result, by the authors 
themselves, is that this common central density may be the 
hallmark of a common formation time, a suggestion which 
appeals to the results of N-body simulations. Navarro et 
al. (1996,1997) found that the central density of halos was 
proportional to the cosmic background density at the very 
earliest stages in their assembly, so halos carry a permanent 
imprint of their formation epoch. 



Confirmation and extension of these ideas has come 
from the combination of N-body simulations with semi- 
analytic models to predict the properties of galaxy popu- 
lations (Cooper et al. 2009; Li et al. 2009; Maccio et al. 
2009; Okamoto & Frenk 2009). All these calculations iden- 
tify the early reionisation of gas and minimum temperature 
for atomic cooling (10 4 K) as the features which lead to com- 
mon central mass in satellites. 

From N-body simulations of three different Milky- Way 
sized host halos, Kravtsov (2009) deduces the central masses 
of subhalos from their accretion epoch and the detailed re- 
sults of higher-resolution simulations. Central mass is found 
to be a power-law in subhalo virial mass, M v . The empirical 
result of constant central mass thus corresponds, in turn, to 
a power-law relationship between virial mass and luminosity. 

These connections between this singular observational 
result and the physics of galaxy formation can be understood 
with particular clarity by investigating the expected distri- 
bution of halos in the central mass - total mass plane. This 
distribution is investigated in f|2| paying particular attention 
to the physical explanation behind each feature and general- 
ising the argument to different measurement radii and host 
halo masses. 

Section[3]takes the argument further to understand how 
the underlying subhalo masses might map on to the lumi- 
nosities of the galaxies contained. Rather than demonstrat- 
ing consistency between these particular results and a par- 
ticular galaxy formation model, the relative influence of in- 
fall, cooling and supernova feedback are explained at a more 
general level. It is then shown how these processes can lead 
to the wide range of luminosities observed. 
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2 THE DISTRIBUTION OF SATELLITE HALO 
MASSES 

2.1 Total halo mass 

The total mass of a satellite halo is defined to be its virial 
mass when it first enters the host halo. After this point, it 
is assumed that its central structure is frozen. Any loss of 
mass would be restricted to the outer shells, not affecting the 
central mass which is of interest here. Such changes (in total 
mass after accretion) are not relevant because the discussion 
concerns the temperature of the halo with regard to cooling, 
which will have mostly taken place while the halo is still 
isolated. (Substructures are not expected to retain their hot 
halos long after accretion). 

The central mass of a given halo is easily calculated by 
assuming the standard dark matter density distribution of 
Navarro et al. (1997) which is parameterised in terms of a 
virial radius, r v , and scale radius, r s : 



f(r) 
/(rv) 



'f{x) = In (l + -) - 



(1) 



where M r is the mass enclosed by radius r and M v = M Tv . 
The ratio of the virial radius to the scale radius is referred 
to as the concentration, c = r v /r s . 

This is a good point to note the definition of virial tem- 
perature, T v (in terms of mean particular mass fimn)'- 



fceTv = 



GM v fj,mH 
2rv 



(2) 



Since this provides a direct relationship between virial radius 
and mass, halos of a given virial temperature will occupy a 
clearly defined locus in the M r — M v plane. 

The subhalo abundance can be predicted, using the 
Extended Press-Schechter theory (Press & Schechter 1974; 
Bond et al. 1991), by building halo merger trees from a 
conditional mass function (Parkinson et al. 2008). At low 
masses, this mass function tends to a simple power-law: 
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as can be seen from the upper panel of Fig. [I] This shows the 
number of satellite subhalos which are predicted to remain 
within a host halo of mass 2.5 x 10 12 Mq, this lying between 
two recent estimate^ for total mass of the Milky Way (Li 
& White 2008; Guo et al. 2009). 

Equation jl| can then be used to find the central mass, 
M r , at some given measurement radius, r. The distribution 
of halos in the M v — M r plane can then be visualised, as 
shown in the main panel of Fig [I] as a series of (solid red) 
contour lines. The projection onto the M r axis is plotted in 
the right panel. 



1 Li & White (2008) found 2.4 X 10 12 M Q by comparing the mea- 
sured recessional velocity (Courteau & van den Bergh 1999) and 
distance (Stanek & Garnavich 1998) of M31 to halos in the Mil- 
lennium simulation (Springel et al. 2005). Guo et al. (2009) found 
2.6 X 10 12 Mq by using the correlation between stellar mass and 
halo mass found in both the Millennuim simulation and the SDSS 
(Li k. White 2009), and applying a value of Af* = 6 X 1O 1O M 
(Flynn et al. 2006). 



2.2 Understanding the distribution of central 
masses 

2.2.1 High mass limit 

The number of halos diminishes rapidly with increasing total 
virial mass fjjp, and there is a corresponding natural upper 
limit for the central mass. However, the exact point at which 
this occurs on the central mass axis is also dependant on the 
choice of measurement radius (300 pc in Figjl]). 

A natural cosmological limit might be expected for the 
quantity of dark matter which can collapse into a region of a 
given size. Such a limit can be estimated by considering the 
characteristic density, p c , of halos, as defined by Navarro et 
al. (1997) in their proposed density profile: 



p(r) 



r/r s (1 + r/r s y 



(4) 



This characteristic density is found to be proportional to the 
mean matter density pm at the very earliest stage^Jin the 
formation of a halo. 



pc ~ 3000pM(fZform) 



3000 af olm f2 M 150M Q kpc" 



(5) 



Choosing, for example, an early formation time with scale 
factor, af orm = 0.1, the relationship |5| gives a characteristic 
density of p c = 10 8 MQkpc -3 . This can be easily translated 
into an enclosed mass by considering the integral of Q in 
the limit r << r a . 



M r « 2nr s r p c or 
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It ought now to be no surprise that the distribution 
of central masses in Fig[l] tails off at M r ~ 1O 8 M0. Also, 
looking at the main panel, it is clear that this high-mass 
limit in M r is not constant in M v . The reason for this is 
clear from the dependence on r s in (JsJ) : this scale radius will 
generally increase with M v , hence the slanted shape of the 
outer contour lines at M v < 10 9 Mq. 

Equations (|5| and (JsJ show that, to acquire a greater 
central mass, halos must have formed even earlier than 
ttform ~ 0.1, or have a much greater scale radius than r s ~ 1 
kpc. The probability of either of these is too low to produce 
a significant population with higher central mass. 



2.2.2 Low-mass limit - baryonic content 

To connect with the observational findings, we must move on 
to consider which subhalos may contain luminous baryons. 
This can be done by reducing the population to include only 
those subhalos whose virial temperature has, at some point 
in their history, exceeded some minimum temperature for 
cooling. 

Such a threshold exists in atomic cooling, at T « 10 4 K. 
Metal line cooling is effective at slightly lower temperatures 
but this will not be a significant process in metal-poor satel- 
lites. Though some early galaxies may form through molec- 
ular hydrogen cooling, the process is inefficient () and we 
assume that its contribution here is negligible. 



2 Specifically, when just 1% of its mass lies in objects exceeding 
10% of the total final mass. 
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Figure 1. The central panel shows theoretical distribution of satellites as a function of their total virial mass, M v , at the point of 
accretion (x-axis) and the mass enclosed within their central 300 pc, M300 (y-axis). Solid (red) lines show all halos. The dashed (blue) 
lines include only the halos whose temperatures (or their progenitors' temperatures; see Fig[2]l have exceeded the threshold for atomic 
cooling. The locus of halos which are at this limit is shown as a thick grey line grey, its width allowing for an appropriately small range 
of halo concentrations, 4 < c < 6 (j |2.1[ l. The upper panel and right panel show the projection of this distribution onto each axis. The 
light blue shading in the right panel highlights the central 50% of the distribution in Af3oo and then, in the top panel, the distribution 
of these same halos in terms of their virial mass. The mass of the host halo at redshift z = is M = 2.5 X 10 12 Mq. 
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Figure 2. The merger tree for two halos in which the final virial 
temperature is insufficient to allow atomic cooling, yet cooling 
was able to occur in their earlier progenitors (red points). Blue 
lines are drawn between the progenitors of each object and a light 
blue point highlights the final object. In both these examples, the 
final cooled baryon fraction, /t,, is extremely low, due to the short 
period of time available for the hot baryons to cool by comparison 
with their free-fall time (j ]3.1| l 



The relevance of re-ionisation to the formation of satel- 
lite galaxies has been highlighted in the articles listed in 
Sjl] After the first stars have formed, background radiation 
is expected to suppress cooling in halos below a given virial 
temperature. Simulations indicate that this limit grows with 
redshift but is about the same order, 10 4 K, as the atomic 
cooling threshold (Okamoto et al. 2008). 

Because the effect of the ionising background is akin to 
that of the limit imposed by atomic cooling, the arguments 
in this article will apply equally well to either of these natural 
thresholds. It is enough to note that the number of "galaxies" 
in the restricted population must always be interpreted as a 
maximum number. 

The locus of T v = 10 4 K is shown in Fig.[l]for an appro- 
priate range of values for the concentration, c. It is noticeable 
that this selected population extends below the cut-off that 
is shaded in Fig[l] The reason for this is as follows. Though 
the halo may have been below the cooling threshold at the 
time of its accretion, one or more of its progenitor halos ex- 
ceeded this limit and may have been able to form a galaxy. 
Two such merger histories from the theoretical subhalo pop- 
ulation are shown in Fig[2] 

As pointed out by many previous authors (Kravtsov 
2009; Li et al. 2009; Maccio et al. 2009) the cooling cut- 
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Figure 3. The distribution in central mass for three different 
measurement radii: 100 pc (bottom row) 600 pc (middle) and 900 
pc (top). In the panels on the left, the contour lines and key are 
the same as for Fig[T] with solid red lines showing contours for the 
entire halo population, and dashed blue lines showing just those 
halos which can contain cooled baryons. The dotted horizontal 
line shows the simple estimate |7]l for the peak in the central 
mass distribution. The panels on the right show the projection of 
this distribution onto the central mass axis. 



off would prevent a swathe of lower-mass halos from cooling 
their baryons. Whilst the total halo population (solid lines in 
the top panel of Fig.[T]) continues to increase at lower masses, 
the reduced, cooled baryon population (dashed lines) tails 
off rapidly at low masses. 

The distribution of central masses for this sub- 
population is dramatically different. When combined with 
the natural limit to the maximum central mass, this leaves 
a remainder distribution sharply peaked around a charac- 
teristic mass, with half the halos occupying just a tenth of 
the total logarithmic range (as shown in the right panel of 
Fig.0. 

This characteristic mass can be understood graphically. 
The halos are "trapped" between two lines, the dotted line 
(r v = 300pc, discussed above) and the curve representing 
the cooling cut-off. This central mass at which these two 
lines intersect gives an approximate value for the resulting 
peak in the distribution: 



A I* 



2rfc B r cu1 
Gfj.mii 



(') 



Evaluating this for the measurement radius of 300 pc gives 
the value: 



2 x 10 7 M P 



(8) 



The physical explanation behind |7| is that the central 
masses of any satellite population will be tightly contained 
within two natural limits. At the high-mass end the limit is 
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Figure 4. The distribution of satellite halos for a range of host 
halo masses. Each of the top four panels shows a single host mass, 
with the predicted number of satellites plotted using the same 
contour lines as Figs.^and[3] Red lines show all halos, and blue 
lines include only those halos whose virial temperature (or their 
progenitors' temperature) has exceeded 10 4 K. The dotted light 
blue line shows the characteristic central mass predicted by equa- 
tion J7j. The lower panel shows the projected distribution on the 
axis for all four host masses and the dotted line again shows 
the characteristic mass. 



set by the maximum dark matter mass which can be enclosed 
by the chosen measurement radius. At the low-mass end it is 
set by the requirement that the halo be large enough to have 
reached a high enough temperature for atomic cooling. This 
latter requirement allows baryons to fall into the centre and 
form a galaxy, which may survive the subsequent traumas 
of accretion into a larger halo and remain as an observable 
satellite. 



2.2.4 Variation of host halo mass 

It is also important to explore the dependance upon host 
halo mass. Changes in this mass will lead to changes in the 
distribution, calculated in Figs, [I] and [3j and it is impor- 
tant to check that the explanation given is not exclusive to 
particular range of halo mass. This is especially interesting 
because the observational constraint^] on the mass of the 
Milky Way's halo are not that tight. 

A range of host halo masses is investigated in Figj4] 
which shows that the argument leading to Q holds equally 
well for all the sub-halo populations. The peak in the distri- 
bution of central mass remains unchanged, the overall pop- 
ulation simply increasing with host halo mass as expected. 



3 THE EXPECTED DISTRIBUTION OF 
LUMINOSITIES 

3.1 Timescales for baryonic infall 

Consider the energy radiation rate per particle, e r and the 
thermal energy per particle, e t h, estimated by making stan- 
dard assumptions: 



eth = -k B T v 



and 



A(T v ,Z)n 



(9) 



The second equation defines the cooling function, A, which 
can be found for given temperature T and metallicity, Z, by 
reference to the published tables of Sutherland & Dopita 
(1993). &b is Botzmann's constant. These two quantities 
provide a cooling timescale for gas at temperature T and 
gas density p: 



t coo i(n,Z,T) = 



eth 



3fc B r 
2n\{T, Z) 



(10) 
(11) 



The value of A is effectively zero below T = 10 4 K but peaks 
at temperatures only slightly higher than this. This means 
that the halo population in which we are particularly inter- 
ested (those which occupy the region just above the cooling 
cut-off) will actually have the shortest cooling times. To al- 
low reference to Fig. |5J we can evaluate (111 at the mean 
density of an isolated halo of mass M v = 3 x 10 8 M Q and 
T v = 2 x 10 4 K, which gives t cool « 4Myr. 

So, if cooling is so rapid, why do these borderline halos 
not cool all their baryons into their central galaxies? The 
time they spend above the cut-off temperature may be short, 
but it comfortably exceeds the cooling time. The explanation 
is that the free-fall time, 



iff 



1 Gyr 



(12) 



2.2.3 The choice of measurement radius 

It is important to widen the discussion and not solely con- 
sider the mass enclosed by one particular radius (so far cho- 
sen to be 300 pc). Equation |7| is based on the central mass 
being tightly constrained by the limits of atomic cooling and 
the maximum mass that can have collapsed into that size of 
region. This argument is as valid for a range of choices for 
enclosure radius as verified by Fig. [3] 



greatly exceeds the cooling time in these cases. So, if a halo 
is fluctuating in temperature as it grows with time, it really 
has to remain above the threshold for a few Gyr if its cen- 
tral galaxy is to contain a large fraction of the halos total 
baryonic mass. 



3 Li & White (2008) give M MW = 2.43 X 10 12 M Q with a 95% 
lower confidence limit of 0.8 X 10 12 Mq. 
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3.2 The importance of outflow 

One process responsible for the ejection of baryons back out 
of the central galaxy is supernovae explosions. This process 
can be modelled by the simple expression for the mass of 
gas which is added to the halo: 



Afout = /3M* 



(13) 



Traditionally (e.g. White & Frenk 1991), and still in many 
recently published models (Bower et al. 2006; De Lucia & 
Blaizot 2007), the "mass- loading" factor f3 is estimated for 
any galaxy on the basis of its maximum rotation speecj^J 
«disk- This is done by assuming that the fraction, esn, of the 
available supernova energy which is carried out of the galaxy 
in outflow is approximately the same for all systems: 



1 ,v 2 esN-EsN „■, 
-M out v disk = Mi, 

2 m SN 



(14) 



Here, Esn is the total energy available from a supernova 
and msN is the mean mass of stars formed per supernova 
(a consequence of the initial mass function). Equations ( 13 1 
and ( |14[ | then yield a physically motivated expression for the 
ratio of outflowing mass to star forming mass: 



{ V hot \ ' 
V Ddisk / 



2 _ 2esN-EsN 

^hot = 

msN 



(15) 



3.2.1 Outflow and final baryon fraction 



Neglecting re-cooling for the moment, the total mass ejected 



at any time is then simply the integral of (13 1, bearing in 



mind the hierarchical formation of the final system, 



Mout(t) = 



^ A (*')Af*i (*')(!*' 



(16) 



E 



Udisk(t) 
Vi(t>) 



Mu(t')dt' 
Mjt) 



(17) 



The subscript, i, runs over all the progenitors of the sys- 
tem which exist at time t' . The variable Vi is a convenient 
abbreviation for the value of «disk f° r progenitor i. 



The last step to (17 1 used some physical constraints 
to pin down the value of the integral. If the smaller, pro- 
genitor systems have lower rotational speeds than the final 
galaxy, the ratio of velocities in square brackets will always 
be greater than one. However, the contribution from each 
stage in assembly is weighted by the fractional contribution 
to the stellar mass, M+dt/M*, so larger systems (with higher 
rotation speeds) will make up a greater part of the integral. 
This will bring the value of the integral close to unity for 
most formation histories. 

Now, in general, the integral ( 16 1 gives the total mass 
ejected over time, not the total mass which stays in the 
hot halo. If the total baryon fraction is close to the cos- 
mological fraction (f2b/^ivi), and the mass of gas retained 
in the galaxy is small compared to the total stellar mass, 
then M out + AL » (il h /Q M ) M v . This leads, with (flTj) , to a 



simple relationship between host halo mass and stellar mass. 

4 This being a suitable proxy for the escape velocity of the sys- 
tem. 



S2m 

"^7 



1 + /3(M* 



(18) 



It is appropriate to write /3 as a function of M* because 
rotational speed is well correlated with stellar mass, both 
observationally (e.g. McGaugh 2005; Kassin et al. 2007) and 
in this modelfj 

The expression ( |18| l can be used to understand the effect 
that the outflow process will have when the distribution of 
satellite halo masses (Figjl]) is translated into a distribution 
of the mass of satellite galaxies. This translation is shown in 
Fig. [5] which plots the relationship ( 18 \ for a range of effi- 
ciencies with which supernova energy is converted to outflow 
(different values of «hot)- 

Included in this figure, to provide some basic reassur- 



ance as to the validity of (18 1, are the equivalent estimates 



derived from an established semi-analytic galaxy formation 
model (Cole et al. 2000), as applied to this very merger tree. 
This solves a full chain of coupled differential equations, one 
of which is (131, that interrelate the various evolving prop- 



erties of a model galaxy populatiorj^] 

This is a good point to redress the first clause of the 
argument supporting ( |18[ ), which asks to neglect re-cooling 
"for the moment". Gas would generally be expected to re- 



cool onto the galaxy on the order of a free-fall time ( 12 1 



However, in this case of satellite systems, the hot gas can be 
stripped off when they are accreted by the host system (Font 
et al. 2008). The question is whether there is usually time, 
before accretion, for a significant fraction of the ejected gas 
to re-cool back on to the galaxy. 

This question can be addressed using the full galaxy for- 
mation model, which includes stripping and recooling. The 
premise, that the former is the dominant effect (Afhot ~ 
Af ou t), is supported by the full model for lower mass satel- 
lites. Larger satellites, on the other hand, can evidently be 
expected to re-cool some of their ejected mass before it is 
stripped. So, the approximation ( 18 1 probably underesti- 



mates the final stellar mass component for larger satellites, 
but still provides a useful analytic guide. 

With its physical aspects qualified, Fig. [5] can now 
be used to understand the effects of feedback on the 
satellite population. Earlier, Fig. [I] identified that 50% 
of galaxy-containing halos lay within a range of just 0.1 
in log(A/3oo/M0). Yet, these same halos were distributed 
across four orders of magnitude in total virial mass. 

When plotted as a function of their baryonic mass, these 
satellites occupy an even greater range due to the effects 
of feedback. This process disperses galaxies from the same 
halo population across seven orders of magnitude in bary- 
onic mass for a high (but reasonable) value of outflow ef- 
ficiency. This provides a plausible explanation for the ob- 
served distribution of the masses and luminosities of Milky 
Way satellites. 

5 Fig[S]uscs i; disk = (A/*/10 10 M Q )°' 22 180kms- 1 , this being a 
good match to the rotational speeds of the model systems with 
which the estimate is being compared 

6 The differential equation relating to star formation (Cole et 
al. 2000, eq. 4.4) might be expected to be crucial here, but this 
mainly controls the cold gas fraction in the galaxies. As galaxies 
with M* < 1O 1O M0 have less than 1% of their mass in cold 
gas, this component can be neglected for the purposes of this 
argument. 
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Figure 5. The variation in baryonic mass resulting from different strengths of supernova feedback. Dashed lines show the expected 
analytic relationship ( |18| between halo mass, M v and the stellar mass contained M*. These different colour show different estimates 
of supernova efficiency, as encapsulated by the parameter v^ ot . when more mass is ejected. The upper panel shows the corresponding 
stellar mass function for each mode. Inset in grey on the lower panel, the mass function for the halos is shown again for easy reference. 
The shaded area shows the halos which occupy the central 50% of the distribution in M300 (see FigQ. The range that these same halos 
occupy on the baryonic mass axis is shown as an error bar for the cases of lowest and highest feedback efficiency. The red arrow highlights 
the increase in the stellar mass of central galaxy (circled) as feedback increases (opposite to its effect on satellite stellar mass). 



3.2.2 Related implications 

Fig[5]also emphasises two concepts about the general effects 
of feedback: 

(i) Outflow from supernovae should not be expected to 
actually reduce the number of satellites, but just reduce 
the fraction of the total baryons from each halo which 
are converted to stars. Higher feedback simply shifts the 
satellite population to lower luminosities, as is clear from 
the baryon mass functions in the upper panel. This shift in 
the satellite stellar mass function reduces the total number 
of baryons which are locked up as stars in the substructure. 



(ii) Though this ejected gas may be stripped and unable to 
return to its original galaxy, it can cool and fall into the cen- 
tre of its new host halo. This has the interesting consequence 
that stronger feedback can actually lead to brighter central 
galaxies, as highlighted in the upper right of the main panel 
of Fig. [5] 
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4 SUMMARY 

We have considered the issue of the present-day central den- 
sity of satellite halos in the context of cold dark matter 
cosmology theory. By generating satellite populations for a 
range of host halo masses, the pattern of their distribution 
in the central mass - total mass plane was understood in 
terms of physical thresholds. 

The central density must be high enough for the satellite 
galaxy to survive the effects of reionisation, and/or for the 
temperature to be high enough to allow atomic cooling. It 
must also be lower than the natural maximum density which 
can arise through gravitational collapse from ACDM initial 
conditions. 

These boundaries restrict the population of satellite 
hosts to a very small range in central mass, but not in total 
mass. The explanation for this effect leads to a simple ana- 
lytic approximation for the central mass |7| given in terms 
of the general values of host halo mass and enclosure radius. 
The merger tree calculations support of this conclusion over 
about an order of magnitude in both factors. 

While half the satellites occupy just a tiny portion of the 
total range in central mass, this same half spans three orders 
of magnitude in pre-infall virial mass. Further analysis then 
helps understand how these virial masses might map onto 
stellar masses, which will explain the range in luminosities 
observed by Strigari et al. (2008). 

By reference to the merger tree formation histories 
(Fig[2| it can be seen that halos may only be able to cool 
effectively for part of their histories. Due to relatively long 
free-fall times, only a fraction of the baryons that are ini- 
tially available in these halos will infall into the galaxies. 
This is one reason why mass-to-light ratios might be high 
for smaller systems. 

A second reason is reheating. In the case of satellites, 
where hot halo gas can be stripped upon accretion, the tra- 
ditional argument (that a set fraction of supernova energy 
is carried out of the galaxy in the reheated gas) is shown to 
lead to a simple relationship between final stellar mass and 
halo mass. 

This makes it clear that even quite modest energy con- 
version will dramatically decrease the fraction of the baryons 
in satellite halos which can be retained in their galaxies. 
For high, but reasonable values of the conversion efficiency, 
the galaxies which occupy a range in central mass of just 
0.1 log 10 (M r /MQ) are distributed across six orders of mag- 
nitude in Mbaryons/MQ. This is in keeping with the wide 
range of luminosities observed for the Milky Way's satel- 
lites. 
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